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The quasi-particle spectra in the metallic rutile and insulating monoclinic phases of VO2 are 
shown to be dominated by local Coulomb interactions. In the rutile phase the small orbital 
polarization among V 3d t2 g states leads to weak static but strong dynamical correlations. In the 
monoclinic phase the large 3d orbital polarization caused by the V-V Peierls distortion gives rise 
to strong static correlations which are shown to be the primary cause of the insulating behavior. 

PACS numbers: 71.20.Be, 71.27.+a. 79.60.Bm 



I. INTRODUCTION 

The metal insulator transition in VO2 has been inten- 
sively studied for a long time. At 340 K the resistiv- 
ity changes by several orders of magnitude^ The high- 
temperature metallic phase has a rutile structure, while 
the low-temperature insulating phase is monoclinic (Mi), 
with a zigzag- type pairing of V atoms along the c axis. 
Both phases are non-magnetic. Although this transition 
is widely regarded as a Mott-Hubbard transitionf 2 *^^ 
the role of the Peierls distortion of the crystal struc- 
ture in the insulating phase has been the topic of intense 
debatei ^'Ti?!?! 1 ?! 11 ! 12 ! 1 ? The discovery of a second insu- 
lating phase (M2)j2 in which only half of the V atoms 
form pairs while the other evenly spaced chains behave 
as magnetic insulators, suggested that both low temper- 
ature phases should be regarded as Mott-Hubbard and 
not as Peierls band insulators^ Although the role of the 
Coulomb interaction in the metal insulator transition of 
VO2 has been studied previouslyji2*ii*i2iiii a consistent 
description of the rutile and monoclinic phases is not yet 
available. 

The aim of the present work is to elucidate the inter- 
play of Coulomb correlations and orbital polarization in 
the quasi-particle spectra of VO2 . Analyzing recent pho- 
toemission datai^i^ we demonstrate that the metallic 
and insulating phases show evidence of strong local inter- 
actions which manifest themselves in distinctly different 
ways because of the different degree of orbital polariza- 
tion in the rutile and monoclinic structures. The key 
aspect of the metallic phase is the small orbital polar- 
ization among V 3d ti g states, implying weak static but 
strong dynamical correlations. Thus, the spectra reveal 
band narrowing and reduced weight of the coherent peak 
near the Fermi level, and an incoherent satellite feature 
associated with the lower Hubbard band. The monoclinic 
phase, in contrast, is characterized by a pronounced t 2g 
orbital polarization induced by the V-V Peierls distor- 
tion. The local Coulomb interaction therefore leads to 
strong static correlations and to the opening of an excita- 
tion gap. The influence of non-diagonal coupling among 
tig orbitals on the size of the gap is investigated within 
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FIG. 1: Measured photoemission spectra for VO2 films on 
Ti02 in the V 3d band region after background subtraction. 
Solid curve: T = 350 K (metallic rutile phase); dashed curve: 
T = 250 K (insulating monoclinic phase); Ef — 0. 14 



the static limit and is found to be small. 

In the following sections we discuss results obtained 
within several theoretical approaches, such as the dy- 
namical mean field theory (DMFT) 1 ^, the local density 
approximation plus Hubbard U (LDA+U) 17 and GWi& 
methods. We argue that none is presently capable of ex- 
plaining all of the observed phenomena in a consistent 
manner. Instead, we focus on the merits and failures of 
these approaches with the aim of highlighting the differ- 
ent roles of Coulomb correlations and orbital polarization 
in the metallic and insulating phases of VO2 . 



II. RESULTS AND DISCUSSION 

Fig. 1 shows photoemission spectra for VO2 films 
(150 A thick) grown on TiO2(001) at binding ener- 
gies corresponding to the V 3d bands (21.2 eV photon 
energy). 14 Emission from O 2p states extends from —2 eV 
to about — 8eV. As a result of compressive stress from 
the substrate 1 ^ the transition temperature is lowered 
from 340 K to 290 K. The high temperature spectrum 
shows a Fermi edge characteristic of metallic behavior 
which is absent in the low temperature spectrum. While 
this trend is in agreement with previous photoemission 
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measurements )22i2i the data in Fig. 1 reveal two spectral 
features in the metallic phase: a peak close to Ep and a 
second one near —1 eV. In the insulating phase, the peak 
near Ep disappears and the feature near — leV becomes 
more intense. Recent VO2 photoemission spectra taken 
at 520 eV photon energy^ agree with the ones shown in 
Fig. 1, except for a considerably greater relative weight of 
the coherent peak near Ep in the metallic phase. These 
spectral changes are consistent with the observation that 
satellite peaks in low photon energy spectra tend to be 
more pronounced as a result of a surface induced en- 
hancement of Coulomb correlations SS&Si 

Before analyzing the photoemission data we discuss the 
single particle properties of VO2 obtained using density 
functional theory. We have carried out full potential lin- 
earized augmented plane wave (LAPW) calculations for 
the rutile and monoclinic structures using the experimen- 
tal lattice parameters and treating exchange correlation 
within the generalized gradient approximation (GGA). 25 
Due to the octahedral crystal field, the states near Ep 
have V 3d ti g character. They are separated by a small 
gap from the empty e g states, and from the O 2p states by 
a gap of about 1.0 eV. The occupancy of the ti g manifold 
is 3d 1 . Our results qualitatively confirm previous LDA 
calculations iSi2ii2iii Fig. 2(a) shows a comparison of the 
V total ti g density of states for the rutile and monoclinic 
phases of VO2 . Although there are differences in detail, 
the overall width of these distributions and the shape of 
the occupied region are similar for both structures. Ev- 
idently the GGA/LDA does not predict the insulating 
nature of the monoclinic phase. Moreover, on the basis 
of these densities one would not expect correlations to 
play very different roles in the two phases. However, if 
we plot the subband contributions to the ti g density, the 
two structures are very different, as shown in Figs. 2(b) 
and (c). Whereas in the rutile phase the ti g bands have 
similar occupation numbers, in the monoclinic phase the 
d x 2_ y 2 band is significantly more occupied than the d xz 
bands. (We adopt the local coordinate system of Ref. 
i.e., x denotes the c axis, while y and z point along the 
diagonals of the a, b plane.) The origin of this orbital po- 
larization is the V-V dimerization along the c axis. The 
d x 2_ y 2 band splits into bonding and anti-bonding compo- 
nents, while the d xz ^ yz bands are pushed upwards due to 
shortening of V-0 distances. Investigation of the energy 
bands shows that, in agreement with earlier workj&£ the 
top of the bonding d x 2_ y 2 bands is separated by a slight 
negative gap from the bottom of the d xz . yz bands. In 
the following subsection we show that different degree of 
orbital polarization in the rutile and monoclinic phases 
has a pronounced effect on the quasi-particle spectra of 
V0 2 . 



A. Metallic Rutile Phase 

Let us first discuss the rutile phase. Comparing the 
tig density of states with the photoemission spectra it is 
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FIG. 2: VO2 3d density of states calculated within LAPW 
method, (a) Total ti g densities for rutile and monoclinic 
phases; (b) and (c) t2 g density of states components for rutile 
and monoclinic phases; Ef = 0. 



plausible to associate the feature close to Ep with emis- 
sion from metallic V 3c? states. The peak near — leV, 
however, lies in the gap between V 3d and O 2p states 
and cannot be understood within the single particle pic- 
ture. To describe the spectra in the rutile phase it is 
clearly necessary to account for dynamical Coulomb cor- 
relations. For the evaluation of the quasi-particle distri- 
butions we use the Dynamical Mean Field Theory com- 
bined with the multiband Quantum Monte Carlo (QMC) 
method. 16 Since hybridization among tig states is weak 
the local self-energy is taken as diagonal in orbital space. 
The tig density of states components shown in Fig. 2(b) 
then serve as input quantities accounting for the sin- 
gle particle properties of the rutile structure. The lo- 
cal Coulomb interaction defining the quantum impurity 
problem is characterized by intra- and inter-orbital ma- 
trix elements U and U' — U — 2 J, where J is the Hund's 
rule exchange integral. According to constrained LDA 
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calculations, U » 4.2 eV and J w 0.8eVii*2& The cal- 
culations are performed at T w 500 K with up to 10 6 
sweeps. The quasi-particle distributions are obtained via 
maximum entropy reconstruction^ 

Fig. 3(a) shows that, in contrast to the single particle 
density of states, the calculated ti g quasi-particle spectra 
for the rutile phase of VO2 exhibit two spectral features, 
a coherent peak near Ep and a lower Hubbard band near 
— 1 eV, in agreement with experiment. The peak near Ep 
accounts for the band narrowing and lifetime broaden- 
ing of the metallic states whereas the Hubbard band is 
associated with incoherent emission. Since the ti g sub- 
bands have comparable single particle distributions their 
quasi-particle spectra reveal similar correlation features. 
Moreover, because of the weak orbital polarization in the 
rutile structure static correlations are negligible. Thus, 
in the metallic phase the spectral weight transfer between 
coherent and incoherent contributions to the spectrum is 
primarily the result of dynamical correlations. 

We point out that, in view of the approximate nature 
of the model underlying the DMFT, quantitative agree- 
ment with photoemission data cannot be expected. On 
the theoretical side, the consideration of purely on-site 
Coulomb interactions and the neglect of the momentum 
variation of the self-energy permit only a qualitative anal- 
ysis of the spectra. In addition, there exists some uncer- 
tainty regarding the precise values of the Coulomb and 
exchange energies. Finally, the DMFT results depend on 
the temperature used in the QMC calculation. The com- 
parison with results obtained for slightly different values 
of U, J and T, however, gives us confidence that in the 
metallic phase the transfer of spectral weight from the 
coherent peak near Ep to the satellite region near —1 eV 
is qualitatively reliable and consistent with analogous dy- 
namical correlation effects in other 3d 1 transition metal 
oxides, such as SrVO^m^^^ As we discuss below, the 
local DMFT treatment predicts the monoclinic phase to 
be also metallic. The metal insulator transition in VO2 
is therefore not achieved simply by lowering the tempera- 
ture in the rutile phase. The lattice transformation from 
the rutile to monoclinic structure must be taken into ac- 
count. Therefore, the DMFT results shown in Fig. 3(a) 
for the rutile structure at T = 500 K can be regarded 
as representative of correlation induced behavior in the 
metallic domain. 

On the experimental side, as pointed out above, it is 
important to recall that photoemission data taken at low 
photon energies represent a superposition of bulk and 
surface contributions. Since correlation effects are ob- 
served to be more enhanced at surfaces ^SiSiSi the rel- 
ative intensity of the satellite peak near — leV in the 
21.2 eV spectra shown in Fig. 1 is considerably more 
pronounced than at high photon energies 15 at which pri- 
marily bulk-like valence states are detected. 

Dynamical effects in the metallic phase may also be 
evaluated within the GW approach^ which treats long 
range Coulomb interactions in the random phase approx- 
imation (RPA) and which has proven rather useful to de- 
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FIG. 3: (a) VO2 ti g quasi-particle spectra for rutile phase 
calculated within DMFT; (b) fo s partial density of states for 
monoclinic phase calculated within LDA+U; (c) monoclinic 
density of states calculated within non-diagonal (solid curve) 
and diagonal (dashed curve) versions of the LDA+U; see text. 
Note the different energy scales. 



scribe excitation spectra of weakly correlated systems. 29 
Because of the neglect of multiple electron-electron and 
hole-hole scattering processes, this scheme fails when lo- 
cal Coulomb interactions are important^ Presumably, 
therefore, for VO2 the GW method does not reproduce 
the lower Hubbard band. The satellite in the rutile phase 
is, of course, also beyond the scope of the static LDA+U 
approach^! The qualitative agreement between the mea- 
sured high-temperature spectra and the theoretical re- 
sults shown in Fig. 2(a) suggests that the DMFT cap- 
tures the key spectral weight rearrangement induced by 
dynamical correlations. 
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B. Insulating Monoclinic Phase 

Turning now to the monoclinic phase we first calcu- 
lated the t2 g quasi-particle spectra within the DMFT. 
Because of the orbital polarization caused by the V-V 
Peierls distortion, correlation effects in the d x 2_ y 2 band 
are now stronger so that the intensity of the lower Hub- 
bard band is enhanced compared to the rutile phase (not 
shown). While this trend agrees with the photoemis- 
sion data, the coherent peak near Ep persists. Thus, the 
single-site DMFT based on a diagonal self-energy does 
not reproduce the insulating nature of the monoclinic 
phased This failure is in striking contrast to the results 
obtained within the LDA+U and GW methodsiiSJ^ 

Let us discuss first the LDA+U approach. In this 
scheme different orbital occupations give rise to orbital 
dependent potential terms which shift different t2 g bands 
in opposite directions. Thus, orbital polarization leads to 
strong static correlation effects. Using the same Coulomb 
parameters as for the rutile phase we find that the lo- 
cal density of states calculated via the LAPW/LDA+U 
exhibits a gap of about 0.7 eV (see Fig. 3(b)), in agree- 
ment with optical data™ and previous LDA+U results— 
The bonding d x 2_ y 2 bands are now completely filled and 
the d xz ^ yz and anti-bonding d x 2_ y 2 bands are shifted up- 
wards. Thus, the LDA orbital polarization in the mono- 
clinic phase is further increased as a result of static cor- 
relations. In turn, this reduction of orbital degeneracy 
enhances the trend towards insulating behavior^ 

We emphasize that the LDA+U does not provide a 
complete description of correlation effects in the insu- 
lating phase. Essentially, it amounts to a fully self- 
consistent, non-local treatment of static screening within 
the dimerized structure. Genuine dynamical effects, such 
as the spectral weight transfer from the coherent to the 
incoherent peak as observed in the metallic phase are ig- 
nored. Neverteless, since static correlations appear to be 
the origin of the excitation gap in VO2, it is instructive to 
inquire which features of the LDA+U contribute to the 
opening of the gap. Of particular interest is the role of 
the non-diagonal occupation matrix n a p which is the key 
input quantity in the orbital dependent perturbation po- 
tential used in the LDA+U. Recent work by Pavarini et 
al*^ on 3d 1 perovskite materials exhibiting non-diagonal 
t2 g orbital coupling caused by octahedral distortions sug- 
gests that this mechanism tends to suppress orbital fluc- 
tuations and to enhance insulating behavior. Fig. 3(c) 
compares the VO2 LDA+U density of states with results 
of an approximate LDA+U treatment in which at each it- 
eration only the diagonal elements retained. The 
size of the gap is seen to be only slightly reduced. In the 
case of VO2, therefore, non-diagonal coupling among t2 g 
orbitals is evidently not the main reason for the existence 
of the gap. 

To understand the gap formation obtained within the 
LDA+U, it is useful to formally express the t2 g self- 
energy matrix as S(w,fc) = E HF (fc) + AE(w,fc), where 
E HF (fc) = H LDA+V (k) - H LDA (k) is real and accounts 



for spectral changes associated with the LDA+U— 
AS(w,fc) is complex and describes purely dynamical ef- 
fects. The important point is that the LDA+U in- 
cludes the full momentum variation of £ HF (fc) within 
the true lattice geometry. Thus, in a site representa- 
tion, static screening processes generate finite inter-site 
elements (each element is a matrix in orbital space) 
even if the bare LDA+U perturbation potential is site- 
diagonal. Thus, S^ F (fc) is more accurate than the static 
limit of the single-site DMFT which neglects the k de- 
pendence and assumes the impurity environment to be 
isotropic, i.e., S^ MFT (o;) ~ §ij. The results shown in 
Fig. 3(b,c) suggest that the proper evaluation of S HF (fc) 
within the Brillouin Zone of the dimerized structure is 
the crucial ingredient to an adequate description of the 
insulating behavior in VO2. The LDA+U amounts to a 
self-consistent treatment of E HF (fc) since the solution of 
the Schrodinger equation imposes no restrictions on how 
the wave functions adjust to the LDA+U potential. 

In order to go beyond static screening and include dy- 
namical correlations in the monoclinic phase a cluster 
extension of the DMFT is most likely required. Such 
an extension is beyond the scope of the present work. 
A cluster DMFT would include the crucial inter-site ele- 
ments S r ^ FT (o;) which arise naturally in a cluster repre- 
sentation of the lattice and of the impurity Green's func- 
tions G{uj) and Go(w), even for a purely on-site Coulomb 
interaction. Preliminary results for VO2 within a clus- 
ter DMFT— show that dynamical screening processes 
beyond the static correlations included in the LDA+U 
cause a broadening of the LDA+U density distribution 
and a shift of the main spectral peak towards the Hub- 
bard bands. Possibly, a multi-site extension of the DMFT 
might also identify the true origin of the metal insulator 
transition in VO2, i.e., whether it is primarily caused by 
the lattice reconstruction or by Coulomb correlations, or 
whether these mechanisms mutually enhance each other. 
All one can say at present is that, given the orbital po- 
larization induced by the lattice transition, Coulomb cor- 
relations have a pronounced effect on the quasi-particle 
spectra. 

As noted above, the GW approach applied to the mon- 
oclinic phase of VO2 also yields an excitation gap of the 
correct magnitude— These calculations utilize a model 
self-energy^ consisting of an approximate short-range 
contribution given by the local exchange correlation po- 
tential, and a correction due to incomplete screening of 
the Coulomb interaction. Essentially, in this simplified 
GW scheme the self-energy correction consists of a "scis- 
sor" operator and additional, non-rigid shifts of energy 
eigenvalues f2S Presumably, the reason why this model 
self-energy yields a gap is that Yj(q,Lu) is non-local and 
non-diagonal in site space, i.e., it includes the important 
static correlations in the dimerized structure in a simi- 
lar fashion as the LDA+U. The approximate nature of 
the GW method, in particular, the neglect of electron- 
electron and hole-hole ladder type interactions, affects 
mainly the remaining dynamical corrections caused by 
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the strong local Coulomb energy. An adequate treatment 
of these corrections would require going beyond the RPA 
and should lead to a more realistic description of the po- 
sition and width of the Hubbard bands. 

An interesting additional feature observed by Okazaki 
et alm^ is the temperature dependence of the photoemis- 
sion spectra below the metal insulator transition. Es- 
sentially towards lower T the excitation gap becomes 
more clearly defined and the main peak near —1 eV gets 
slightly sharper. These changes, however, are very small 
on the scale of the main discrepancy still existing between 
the LDA+U or GW results and the experimental data. 
At present it is not clear whether cluster DMFT calcula- 
tions in the monoclinic phase as a function of tempera- 
ture will be able to explain the observed trend or whether 
an explicit treatment of electron-phonon coupling is re- 
quired. 

III. SUMMARY 

The metal insulator transition in VO2 appears to be 
remarkably complex and its origins are not yet fully un- 
derstood. In the present work we focussed on the im- 
portant role of two aspects, local Coulomb correlations 



and orbital polarization, in the low and high tempera- 
ture photoemission spectra of VO2. Whereas the metal- 
lic phase exhibits weak static and strong dynamical cor- 
relations, the monoclinic phase is dominated by static 
Coulomb correlations. Accordingly, the rutile spectra re- 
veal a double-peak structure where the feature close to 
Ep is identified with metallic V 3d states and the peak 
near — leV with the lower Hubbard band. Since the ti g 
states in the metallic phase are roughly equally occu- 
pied orbital polarization is negligible. The fundamental 
difference in the monoclinic phase is the large orbital po- 
larization induced by the symmetry breaking due to V-V 
dimerization. The preferential occupation of the d x 2_ y 2 
bonding states implies strong static correlations which 
are the main origin of the excitation gap. It would be 
of great interest to study the additional dynamical cor- 
relation effects in this phase within an extension of the 
single-site DMFT approach. 
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